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Abstract. Star-disc coupling is considered in numerical models where the stellar field is not an imposed perfect dipole, but 
instead a more irregular self-adjusting dynamo-generated field. Using axisymmetric simulations of the hydromagnetic mean- 
field equations, it is shown that the resulting stellar field configuration is more complex, but significantly better suited for 
driving a stellar wind. In agreement with recent findings by a number of people, star-disc coupling is less efficient in braking 
the star than previously thought. Moreover, stellar wind braking becomes equally important. In contrast to a perfect stellar 
dipole field, dynamo-generated stellar fields favor field-aligned accretion with considerably higher velocity at low latitudes, 
where the field is weaker and originating in the disc. Accretion is no longer nearly periodic (as it is in the case of a stellar 
dipole), but it is more irregular and episodic. 
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1. Introduction 



Over the past decade, numerical simulations have been used 
to address the problem of star-disc coupling by a stellar dipole 
magnetosphere (Hayashi et al. 1996, Hirose et al. 1997, Miller 
& Stone 1997, Goodson & Winglee 1999, von Rekowski & 
Brandenburg 2004, hereafter referred to as vRB04, Romano va 
et al. 2004). One of the perhaps most striking differences 
compared with the standard pictures of star-disc coupling 
(Konigl 1991, Cameron & Campbell 1993, Shu et al. 1994, 
Hartmann et al. 1994) is the realization that the stellar mag- 
netospheric field may not permeate the surrounding disc over 
broad enough a range in radius to ensure sufficient coupling 
between the two. Instead, differential rotation between disc 
and star leads to magnetic helicity injection that in turn in- 
flates the closed magnetic loops which then open up close to 
the inner disc edge (Lovelace et al. 1995, Agapitou & Pa- 
paloizou 2000, Uzdensky et al. 2002). This led Matt & Pu- 
dritz (2004, 2005) to the conclusion that the star cannot be 
braked by the magnetic interaction with the disc, as was as- 
sumed so far. Instead, stellar winds may be the more impor- 
tant agent causing stellar spin-down. 
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One of the important new features used in the simula- 
tions of vRB04 is the fact that the disc itself can produce a 
magnetic field. This feature was introduced by von Rekowski 
et al. (2003) to study the possibility of collimated outflows 
in the absence of an ambient magnetic field. This idea seems 
now quite reasonable also from an observational point of view. 
Menard & Duchene (2004) report observations of the Taurus- 
Auriga molecular cloud which is a nearby star-forming re- 
gion of low-mass stars. They find that jets and outflows from 
protostellar star-disc systems in this cloud are not always 
aligned with the local magnetic field of the cloud where they 
are born. This suggests that jets are not likely to be driven 
by the cloud's own field. However, theoretical models have 
shown that magnetic fields are important for jet acceleration 
and collimation. Further, there is now also observational ev- 
idence that jets and outflows are essentially hydromagnetic 
in nature. Observations of CO outflows exclude purely radia- 
tive and thermal driving of jets and outflows from low-mass 
protostellar star-disc systems (Pudritz 2004). Thermal pres- 
sure might still be dynamically important to lift up the winds. 
Moreover, there are observational indications of warm wind 
regions with temperatures of a few times 10 4 K (Gomez de 
Castro & Verdugo 2003). These winds cannot be due to stel- 
lar radiation because the star is too cool, suggesting that the 
warm winds are magnetohydrodynamically (MHD) driven. 
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So it is quite plausible that a disc dynamo might be necessary 
for the launching of jets and outflows. 

Also the protostar's own magnetic field could contribute 
to the launching of jets and outflows. From observations it 
is not clear whether the jets and outflows originate from the 
central star or from the surface of the inner part of the cir- 
cumstellar disc. In the present paper we extend the simula- 
tions of von Rekowski et al. (2003) and vRB04 by allowing 
the stellar magnetic field to be self-generated by its own dy- 
namo and to react to its environment. One of the immediate 
consequences of a stellar dynamo is that the stellar field is 
no longer an ideal dipole, but it has a range of higher multi- 
poles. Indeed, observations of solar- type low-mass protostars 
surrounded by a disc (classical T Tauri stars, hereafter CTTS) 
suggest kilogauss fields that are rather irregular and in general 
nonaxisymmetric (Johns-Krull et al. 1999a,b, Johns-Krull & 
Valenti 2001). Valenti & Johns-Krull (2004) made measure- 
ments of the magnetic field strength and geometry at the sur- 
face of CTTS. They find a mean surface field strength of 
about 1 kG to 3 kG and a top surface field strength of as large 
as about 6kG. Their measurements indicate that the stan- 
dard magnetospheric accretion model (where a stellar dipole 
is threading the disc leading to a funnel-shaped polar accre- 
tion flow along the magnetospheric field lines) seems far too 
simple and idealized. Their observations show that not a fixed 
fraction of the stellar surface is coupled to the disc by a mag- 
netic field, and that the magnetic star-disc coupling is also 
not periodic in time, both opposed to the standard magneto- 
spheric accretion models. Accretion and the stellar magne- 
tosphere are rather irregular and nonaxisymmetric, and they 
rule out a global dipolar field at the stellar surface of at least 
some CTTS. Although CTTS seem to have complex mag- 
netic fields that are not simple dipoles, the stellar fields still 
seem to have a very large scale component in and around the 
surface of at least a few hundred gauss. 

Observations (e.g. with the Hubble Space Telescope) show 
that outflows from protostellar star-disc systems can have 
different shapes. They can be weakly collimated gas bub- 
bles or well-collimated jets extending far into space. Outflow 
speeds range between 50 and 500 km/s. Interestingly, high- 
resolution observations with the VLT show that also the ac- 
cretion flow velocity is rather high, of the order of 100 km/s 
(Stempels & Piskunov 2002). Protostellar systems are very 
dynamic, with changes in the disc brightness, and knots and 
gaps in the jet. Stempels & Piskunov (2002) find that the time 
variability related to accretion and outflow processes occurs 
on time scales ranging from years to minutes. 

In this paper, we investigate the effects of a stellar dy- 
namo on star-disc coupling and outflows. The accretion flow 
geometry is not known from observations other than it seems 
to be complex and highly time-dependent. Our stellar dy- 
namo models suggest irregular and loose magnetic star-disc 
coupling by the disc dynamo-generated field at low latitudes. 
Accretion is episodic but fast, and along field lines in those 
equatorial regions, where the field is very weak. 

We further suggest that the observed jets and outflows 
could originate from the star and/or from the inner disc. In 
the presence of stellar and disc dynamos, stellar and disc wind 



mass loss rates are comparable, and stellar wind speeds even 
exceed disc wind speeds. The fast stellar wind significantly 
contributes in braking the star. 

2. The models 

A detailed description of the models without stellar dynamo 
and of the models with an anchored stellar dipole magneto- 
sphere can be found in Sect. 2 of von Rekowski et al. (2003) 
and vRB04, respectively. We begin the description of the present 
models by reviewing the basic setup of the earlier models. 
The implementation of the stellar dynamo is described in 
Sect.lO 

2.1. Basic equations 

We solve the set of MHD equations using cylindrical polar 
coordinates (tu, ip, z). In all our models, we solve the conti- 
nuity equation, 

|f + V-(^)=9f c + < ar , (1) 
the Navier-Stokes equation, 

W = ~\ VP ~ V$ + ~gt F + (UKep ~ u)q e SC } ' (2) 

and the mean-field induction equation, 

dA 

— — = wxB + aB — rinoJ. (3) 
at 

Here, u is the velocity field, g is the gas density, p is the gas 
pressure, is the gravitational potential, t is time, D/Dt = 
d/dt+u-V is the advective derivative, F = Jx V-r v i SC 
is the sum of the Lorentz and viscous forces, J — V x B / j^o 
is the current density due to the mean magnetic field B, {1q 
is the magnetic permeability, and r v i SC is the viscous stress 
tensor. To ensure that B is solenoidal, we solve the induction 
equation in terms of the magnetic vector potential A, where 
B = VxA. 

In the disc, we assume a turbulent (Shakura-Sunyaev) 
kinematic viscosity, v t — assc s ^o, where ass is the Shakura- 
Sunyaev coefficient (less than unity), c s = (jp/g) 1 ^ 2 is the 
sound speed, 7 = c p /c v is the ratio of the specific heat at 
constant pressure, c p , and the specific heat at constant vol- 
ume, c v , and zo is the disc half-thickness. 

The mass in the star and the disc that is explicitly ac- 
counted for in the simulations, cannot include the contribu- 
tions from the most dense inner parts. For this reason we 
have, as described in vRB04, modeled this by a self-regulatory 
mass source in the disc and a self-regulatory mass sink in the 
star, q^ lsc and ^ tar , respectively. Matter is injected with Ke- 
plerian speed i/Kep in the disc and removed from the central 
parts of the star without affecting its angular velocity. 

Star, disc, and corona are represented as regions with three 
different poly tropes, but the same polytropic index. The method 
has been described in detail by von Rekowski et al. (2003) 
and vRB04, and is summarized in Appendix | Appendix A:| 
Axisymmetry is assumed throughout this paper. 
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2.2. Disc dynamo and stellar dynamo 

To date, almost all models of the formation and collimation 
of winds and jets from protostellar accretion discs rely on 
externally imposed poloidal magnetic fields and ignore any 
dynamo-generated poloidal field produced in the disc or in 
the star. This is not the case in the models developed in our 
earlier papers, which also form the basis of the models used in 
the present paper. In von Rekowski et al. (2003) and vRB04, 
we have studied outflows and accretion in connection with 
magnetic fields that are produced by a disc dynamo, resolving 
the disc and the star at the same time as well as assuming 
a non-ideal corona. The seed magnetic field is large scale, 
poloidal, of mixed parity, weak and confined to the disc. The 
mixed-parity choice allows for the symmetry to be violated 
also in the (statistically) saturated magnetic fields. 

Our present models are similar to those of von Rekowski 
et al. (2003) and vRB04, except for the modeling of the stellar 
magnetic field. In von Rekowski et al. (2003), the magnetic 
field in the star-disc system is solely generated and main- 
tained by the disc dynamo so that the stellar field is produced 
entirely by advection of the disc field. In vRB04, we also 
model - in addition to the disc dynamo - a stellar magne- 
tosphere that is initially a dipole threading the disc but can 
freely evolve with time outside the anchoring region starting 
at about 1 .5 stellar radii. In the following we usually refer to 
this type of model as "stellar dipole" model (or similar) al- 
though it is important to note that the dipole is imposed only 
in the anchoring region (which includes the star and extends 
up to 1.5 stellar radii). 

In this paper we assume in our main models that both the 
disc field and the stellar field be generated by disc and stellar 
dynamos, respectively. These dynamos are considered to be 
the only mechanism for the generation and maintenance of 
the entire magnetic field in our system. We study the interac- 
tion of both stellar and disc fields, the resulting field config- 
uration as well as the associated outflow and accretion pro- 
cesses. We assume that the magnetic fields in the disc and in 
the star be generated by standard a 2 Vt dynamos (e.g., Krause 
& Radler 1980), where a is the mean-field a effect and Vi is 
the angular velocity of the plasma. This implies an extra elec- 
tromotive force, aB, in the induction equation for the mean 
magnetic field, B, that is restricted to the disc and to the star. 
As usual, the a effect is antisymmetric about the midplane. 
We take 

Ql — -U (t-) 

zo 1 + «i/( C 2MaL c ) r, 1 + ^/(c^MaLr) ' 
where va is the Alfven speed based on the total magnetic field 
(i.e. va = \B\/y/fiog) 9 £disc and £ sta r are time-independent 
profiles specifying the shapes of the disc and the star, respec- 
tively (see Appendix | Appendix A:| ), r* is the stellar radius, 
and QfQ lsc and aQ tar are a effect coefficients that control the 
intensity of dynamo action, whereas Madi SC and Ma star are 
parameters that control the strength of the saturated magnetic 
field. 

In the disc, we choose the a effect to be negative in the 
upper half (i.e. aQ lsc < 0), consistent with results from sim- 
ulations of accretion disc turbulence driven by the magneto- 
rotational instability (Brandenburg et al. 1995). The resulting 



dynamo-generated magnetic field symmetry is roughly dipo- 
lar because the accretion disc is embedded in a (not perfectly) 
conducting corona rather than a vacuum (see vRB04 for more 
references). In the star, we choose the a effect to be positive 
in the upper hemisphere (i.e. aQ tar > 0), as is expected if tur- 
bulence in the star is due to convection. Further, we include 
a quenching which leads to the stellar and disc dynamos sat- 
urating at a level close to equipartition between magnetic and 
thermal energies, depending on Madi SC and Ma star . 

Note that by letting the a effect operate in the whole 
star, we implicitly assume our protostar to be fully convective 
which can indeed be the case for classical T Tauri stars. Cor- 
responding three-dimensional turbulent MHD simulations by 
Dobler et al. (2005) have shown that outside a fully convec- 
tive star, the stellar dynamo-generated magnetic field is dom- 
inated by the large scale field that has a strong poloidal com- 
ponent with a significant dipole moment. They also discuss 
the effect of having no stably stratified overshoot layer, that 
is normally thought important for producing large scale fields 
(Durney et al. 1993, Hawley et al. 2000). However, buoyant 
magnetic flux tubes can remain within the convective star due 
to the effect of turbulent pumping, as is seen in simulations of 
stratified convective dynamos (Nordlund et al. 1992, Tobias 
etal. 1998). 

We recall that the star is rotating differentially with re- 
spect to both radius and latitude, as a consequence of the ini- 
tial setup and determined mainly by the magnetic feed-back 
as time evolves. The star belongs entirely to the computa- 
tional domain. In our stellar dynamo models, no boundary 
conditions are imposed on the velocity or magnetic field at the 
stellar surface. The resulting dynamo effect due to the differ- 
ential rotation (ft effect) is difficult to quantify in and around 
the star. The magnetic diffusivity rj is finite in the whole do- 
main, and dependent on position. 

On the outer boundaries of the computational domain, on 
w = tu max and z = ±£ max , we impose outflow conditions 
(see also Sect. l2.3L On the axis (w = 0) we impose regularity 
conditions. 

2.3. Normalization and model parameters 

In all our simulations we use dimensionless variables. Our 
models can be rescaled easily, and thus they can be applied 
to a range of different astrophysical objects. Here, we con- 
sider values for our normalization parameters that are typi- 
cal of a protostellar star-disc system. As in vRB04, we scale 
the sound speed with a typical coronal sound speed of c s o = 
10 2 kms _1 , which corresponds to a temperature of To ~ 
4 x 10 5 K, and we scale the surface density with a typical 
disc surface density of So = 1 gem -2 . Further, we assume 
M* = 1 Mq (where M* and M are the stellar and solar 
mass, respectively), a mean specific weight of /i = 0.6, and 
7 = 5/3. Fixing our normalization parameters (c s o or To, 
and So, M*, /i and 7) defines the units for all quantities. [In 
the following, G w 6.67 x 10 -8 cm 3 g _1 s -2 is the grav- 
itational constant and 1Z w 8.3 x 10 7 cm 2 s _2 K _1 is the 
universal gas constant. With M w 2 x 10 33 g, it follows 
that GM* = GM Q « 1.3 x 10 26 cm 3 s" 2 .] The result- 
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Table 1. Varying model parameters of the runs discussed in this paper. Note that Model N of vRB04 is different from Model N 
of the present paper in that in vRB04 the gap between star and disc is much smaller and the magnetic diffusivity is as large as 
in Model S there. In Model N, there is no stellar dynamo or anchored magnetosphere but there is a disc dynamo. Model S is 
very similar to Model S of vRB04 (with an anchored stellar magnetosphere). The difference is that in this paper the angular 
velocity profile in the star is not fixed, whereas in vRB04 it was fixed to rather small values. (The asterisk in Model S 
indicates that the magnetosphere is anchored up to a radius of about 0.275, whereas the stellar radius is located at 0.15.) The 
main differences between the models of the runs discussed in this paper are marked in bold in this table. Model Ref is the 
reference model of this paper, with a stellar dynamo as well as a disc dynamo. Model SmS has a smaller star, and Model SDd 
has a stronger disc dynamo, both compared to the reference model. Model NDd has no disc dynamo, and Model ND has no 
disc at all. In order to produce a setup without a disc, the following model parameters are changed (marked with ** in the 
table): (3^ sc = 1 for the specific entropy, aQ lsc = and no seed magnetic field in the original disc for the dynamo, ass = 
and ?7disc,turb = for the turbulent diffusivities, and no mass source. 



Model 


Description 


r* 


/^disc 


„ disc 


Ma d isc 


a s tar 


Ma sta r 




N 


No stellar dynamo, no magnetosphere, but disc dynamo 


.15 


.005 


-.15 


1 










SmS 


Smaller Star compared to reference model 


.15 


.005 


-.15 


1 


+1.5 


10 





Ref 


Reference model with stellar and disc dynamos 


.25 


.005 


-.15 


1 


+1.5 


5 





S 


Strong anchored magnetosphere and disc dynamo 


.15- .275* 


.005 


-.1 


3 







25 


SDd 


Stronger Disc dynamo compared to Model Ref 


.3 


.005 


-.3 


2 


+1 


5 





NDd 


No Disc dynamo but stellar dynamo 


.3 


.005 







+1 


5 





ND 


No Disc but stellar dynamo 


.25 


1** 


o** 




+1.5 


5 






ing units for our dimensionless quantities are given in Ap- 
pendix |Appendix B:| 

The axisymmetric computations have been carried out in 
the domain {w,z) G [0, 2] x [—1, +1], with the mesh sizes 
bw = Sz = 0.01. In our units, the domain then extends to 
±0.1 AU in z and 0.2 AU in w. [In Model ND, the domain 
is larger in z, with —2 < z < +2; the mesh sizes are the 
same.] We impose regularity conditions on the axis (w = 0) 
and do not allow for inflows at the outer boundaries of the 
computational domain, on w = tu max and z — ±£ max . The 
equations are solved also on the boundaries; the derivatives 
are calculated with single- sided schemes of the same order 
(sixth order) as in the domain. 

We choose a set of values for our model parameters such 
that the resulting dimensions of our star-disc system and the 
resulting initial profiles of the physical quantities are close 
to those for a standard accretion disc around a protostellar 
object. Our (main) model parameters are the stellar radius, 
disc inner radius and disc height describing the geometry of 
the star-disc system, the entropy contrast parameters /5di S c 
and /? s tar> me (turbulent) diffusivity coefficients, and the a 
effect parameters for the disc and for the star. 

We vary the stellar radius between r* = 0.15 and r* = 
0.3, corresponding to 3 to 6 solar radii. The disc geometry 
remains unchanged, with the disc inner radius vo- in = 0.6 
(which is therefore located at 4 to 2 stellar radii) and a disc 
semi-thickness of zq = 0.15 (i.e. 1 to 0.5 stellar radii). The 
disc outer radius tu out is close to the outer domain boundary. 

In our piecewise polytropic model, we need to choose 
values for the polytrope parameters in the disc and in the 
star (see Appendix | Appendix A:| ». Our choice for the entropy 
contrast between disc and corona, /^disc = 0.005, leads to an 
initial disc temperature ranging between 9000 K in the inner 
part and 900 K in the outer part. Real protostellar discs have 
typical temperatures of about a few thousand Kelvin (e.g., 



Papaloizou & Terquem 1999). As turns out from our simu- 
lations, the disc temperature increases by less than a factor 
of 2 with time, except for the inner disc edge where the in- 
crease is much higher. The low disc temperature corresponds 
to a relatively high disc density. In the inner part of the disc, 
it is about 3 x 10 -10 g cm -3 initially, but increases to about 
10 -9 g cm -3 at later stages. In the outer part it is a few times 
10 -11 g cm -3 . For the star-corona entropy contrast we choose 
always f3 star = 0.02. 

The Shakura-Sunyaev coefficient of the turbulent disc kine- 
matic viscosity, v t = assc s £o£disc(^ z), is ass = 0.004 in 
all models with disc, whereas the turbulent disc magnetic dif- 
fusivity rjt = 77disc£disc(^ z) is also dependent on position 
but varies in the different models. In the stellar dynamo mod- 
els, agg = 7?disc/ (c s ^o) is typically a few times higher than 
(n) 

ass, but agg can also be as much as an order of magnitude 
higher than ass- Values of the effective magnetic diffusivity 
in the disc range from around 0.0003 in the outer parts to top 
values of around 0.003 toward the inner disc edge so that in 
the disc, the diffusion time tdiff = z o I ^disc ranges between 
75 and 7.5 in the stellar dynamo models. In the stellar dipole 
models, the effective magnetic diffusivity in the midplane is 
a few times 10 -3 , leading to tdm ~ 5 ... 10. 

The dimensionless dynamo numbers related to the a ef- 
fect in the disc and in the star, 1Z^ 1SC = aQ lsc zo/^disc and 
ft s a tar = ag tar r* faster, are7^ isc w -75 . . .-7.5 and7^ s a tar « 
+100 ... + 200 in Model Ref, and 7^ isc « -3 ... - 7 in 
Model S. Note that ag isc « (1 ... 3) c s in Model Ref and 
Model S, and ag tar ~ 2 c s (near the stellar surface) in Model Ref. 
Since CTTS can be fully convective, the a effect would op- 
erate over a comparatively large volume so that in these pro- 
tostars the dynamo may be operating in a highly supercritical 
regime. 

With our initial setup as the hydrostatic equilibrium of a 
piecewise polytropic model with smoothed gravitational po- 
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Fig. 1. The emergence of a fast stellar wind in the presence of a stellar dynamo (right hand panel) compared to the case 
without stellar dynamo (left hand panel) where the stellar wind is very weak. Black arrows: poloidal velocity vectors; solid 
lines: poloidal magnetic field lines (white means right-handed and red/gray means left-handed); colors/gray shades: tempera- 
ture T [log 10 h = (-2, -1, 0, 1) corresponds to T « (3x 10 3 , 3x 10 4 , 3x 10 5 , 3x 10 6 ) K]; black dashed line: poloidal Alfven 
surface. Left: Model N (where c^ tar = 0) at time t « 940 days (steady state). Right: Model SmS (where c^ tar = +1.5) at 
time t « 1644 days (steady state for the disc dynamo). 



tential, the initial stellar surface angular velocity depends on 
the chosen stellar radius. In our models, this angular velocity 
varies at the equator between 1 (for r* = 0.15) and 0.1 (for 
r* = 0.3), corresponding to rotation periods of 10 days to 
100 days. In the case of a rotation period of around 10 days, 
the corotation radius is vo co « 1 in nondimensional form. The 
inner disc radius (which is fixed in our models) is rotating 
with roughly Keplerian angular velocity, resulting in a rota- 
tion period of the inner disc edge of around 4.5 days in all our 
models. The Keplerian speed of the inner disc edge is about 
130km/s, and the free-fall velocity is = y/2GM*/r w 
183 km/s ... 283 km/s in the gap, between the inner disc ra- 
dius and the stellar radius of Model Ref . The Keplerian or- 
bital period of the disc radius close to our outer boundary is 
about 25 days. For an accretion velocity of 10 km/s, the ad- 
vection time through the disc is t accr = (tu ou t — ^in) /^accr ~ 
20 days. We ran the simulations discussed in this paper over 
many rotation periods, between 165 and 1650 days. 

For the stellar and accretion disc dynamos to work, we 
have weak seed magnetic fields in both the disc and the star. 
There is no externally imposed magnetic field in the system in 
all our new models presented in this paper. The initial hydro- 
static solution is an unstable equilibrium because of the ver- 
tical shear between the disc, star and corona (Urpin & Bran- 
denburg 1998). In addition, angular momentum transfer by 
viscous and magnetic stresses - the latter from the stellar and 
disc dynamos - drives the solution immediately away from 
the initial state. The initial state should certainly not be re- 
garded as an approximation to the final time-dependent state. 



In fact, all the results presented in the next section are rep- 
resentative of a statistically steady state, reached after initial 
transients have died out. 

3. Results 

A summary of the varying input parameters of the models 
discussed in this paper is given in Table [l] They will be ex- 
plained in more detail below as the models are being moti- 
vated. We begin with Model N which has no stellar dynamo 
and no anchored magneto sphere, but it has a disc dynamo. 

3.1. Stellar dynamo versus no own stellar field 

In Fig.[I]we show the resulting flow and field lines of Model N 
and compare with those of Model SmS where there is a stel- 
lar dynamo. In both models, the stellar radius (about 3 solar 
radii) and the inner disc edge (about 12 solar radii) are the 
same. Obviously, this comparison is somewhat artificial, be- 
cause in a model without a stellar dynamo the disc would in 
reality be truncated at a smaller radius. Here we have chosen 
the same gap, because we want to see the effect of a stel- 
lar dynamo while keeping the other parameters the same. We 
will discuss a smaller gap in the case without a star's own 
magnetic field at the end of this subsection. 

Comparing the left and right hand panels of Fig.[l] it be- 
comes clear that a stellar dynamo affects the dynamics of a 
star-disc system significantly. In Model SmS, the stellar dy- 
namo generates a magnetosphere that is mostly of roughly 
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Fig. 2. The stellar dynamo-generated magnetic field switches between dipolar (left hand panel) and quadrupolar (right hand 
panel) symmetry. Shown is the same as in Fig.[T] but for Model Ref where the stellar radius is larger. Left, at an earlier time 
(t « 1171 days) when the stellar field configuration turns out to be close to dipolar. Right: at a later time (t w 1183 days) 
when the stellar field configuration turns out to be close to quadrupolar. 



dipolar symmetry and has a magnetic field strength of £? s tar,mc 
300 G. As a result, the star is rotating fast with an equato- 
rial surface rotation period of P ro t,surf ~ 2 days, and a fast 
stellar wind is launched with a terminal speed of about 500 
km/s. This fast wind from the stellar dynamo is supersonic 
and super- Alfvenic and tends to collimate at low heights; it 
is not present in the absence of a stellar dynamo (Model N) 
where the stellar wind reaches 20 km/s along the rotation 
axis and 5 km/s further away from the axis, and where the 
star rotates about 50 times more slowly than in Model SmS. 
In Model N, the stellar wind is subsonic and sub- Alfvenic. 
When there is no dynamo in the star, the stellar field results 
entirely from advection of the dynamo-generated disc field, 
which is only a few tens of gauss strong. In both cases, the 
compressed magnetic field at the inner disc edge acts like a 
shield and prevents any pronounced accretion. 

The disc wind is robust and its overall structure is largely 
independent of the origin of the stellar magnetic field. The 
disc wind terminal speed is about twice the Keplerian speed 
of the inner disc edge. However, the stellar dynamo and in- 
duced stellar wind distort the disc dynamo and disc wind 
in places, and the terminal speed of about 350 km/s in the 
undistorted inner disc wind of Model N drops to about 250 
km/s in Model SmS. Therefore, the total wind mass loss rate 
^wind is dominated by the dense high-speed stellar wind in 
the model with stellar dynamo, where M w ind ~ (2 ... 6) x 
1O _T M yr _1 , and by the less dense but fast disc wind in the 
model without stellar dynamo, where M w ind ~ (1.5 ... 2) x 
10- 7 M o yr- 1 . 

We now want to discuss the effect of a smaller gap in 
the case without the star's own magnetic field. To this end, 



we compare Model N of this paper and Model N of vRB04. 
The smaller gap in Model N of vRB04 has significant con- 
sequences. The smaller gap leads to an initially much faster 
star, with a rotation period of a few days. Further, due to the 
small gap, the advected disc magnetic field is now penetrating 
the star and it thus creates a stellar field strength of more than 
100 G. This fast magnetized star is the reason why the stellar 
wind velocity and mass loss rate are higher in Model N of 
vRB04. Also the inner disc edge - where the strongest disc 
wind originates from - is rotating faster as it is further inward. 
Therefore, also the disc wind velocity and mass loss rate are 
higher in vRB04. Finally, the mass accretion rate is signifi- 
cantly larger in Model N of vRB04. This is because there, the 
build-up of a strong magnetic field advected from the disc oc- 
curs only in the star, whereas in Model N of this paper there is 
a field in the gap, where it is built up to about 250 G and acts 
as a barrier for the disc matter. (The field strength in Model N 
of this paper decreases then to a few tens of gauss in the star.) 

3.2. Dipolar and quadrupolar stellar dynamos 

As our reference model, we choose a model where the stel- 
lar radius is increased to about 5 solar radii. The inner disc 
radius remains at the same position (about 12 solar radii) so 
that it is now located at 2.4 stellar radii, whereas in the mod- 
els of Fig.[l]it is located at 4 stellar radii. By increasing the 
stellar radius without changing the disc location or geometry, 
this results in a more realistic aspect ratio between the stellar 
radius and the height of the inner disc edge and, furthermore, 
the gap and dynamo-free region are decreased so that star- 
disc interaction is facilitated. In Fig. [2] we show the resulting 
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Fig. 3. Space-time diagram for Model Ref in a color/gray 
scale representation of the radial magnetic field component. 
No clear latitudinal migration pattern of the stellar surface 
field is seen, although some high-latitude structures appear 
to show traces of poleward migration. Bright colors or light 
shades indicate positive values; dark colors or dark shades in- 
dicate negative values. The maximum strength of B r is about 
180 G. The vertical lines at times 762.5 and 770 mark the cor- 
responding times t « 1171 days and 1183 days, respectively, 
i.e. the times of Figs|2|and|4l 



axisymmetric spherical a 2 dynamo (Meinel & Brandenburg 
1990). Tavakol et al. (1995) studied the structural stability of 
axisymmetric a£l dynamos in spheres subject to the increase 
of TZ a . For negative radial rotational shear (as is also the case 
in some of our models; see Sect. [321 below), their stable so- 
lution switches from antisymmetric to mixed parity to oscil- 
lating symmetric to symmetric. Our simulations have been 
carried out with a 2 Q dynamos with large 7^^ tar in the star. 
The resulting irregular behavior of the magnetic field is in ac- 
cordance with observations of CTTS, which show that CTTS 
have highly time-dependent fields suggesting the existence of 
irregular magnetic cycles (cf. Johns-Krull et al. 1999a,b). 

The maximum stellar field strength varies also with time, 
ranging between 250 G and 750 G inside the star. The rota- 
tion period of the stellar surface is oscillating around a few 
days around the equator, and the stellar wind reaches termi- 
nal speeds between 300 km/s and 400 km/s; it becomes su- 
personic and super- Alfvenic very quickly. The disc field is 
less than 50 G strong in the bulk of our disc, but it can grow 
to over 100 G around the inner disc edge. The disc wind is 
little affected. It is fastest within a conical shell originating 
from the inner disc edge and reaches a maximum terminal 
speed of up to 250 km/s. The combined stellar and disc wind 
^ (2 ... 4) x 1Q- 7 M yr-\ while 



magnetic field and flow patterns. The accretion disc dynamo- 
generated magnetic field is - as usual - of dipolar symmetry 
with respect to the disc midplane (within the disc as well as 
above and below it), with < in the conical shell above 
the disc and > in the conical shell below the disc. 
The stellar dynamo generates a magnetosphere that switches 
between 'dipolar' (Fig. left hand panel) and 'quadrupo- 
lar' (Fig. right hand panel) symmetries (see also Fig. [3]). 
Here, the quotation marks on 'dipolar' and 'quadrupolar' in- 
dicate that, within the star, both poloidal and toroidal fields 
have dipolar/quadrupolar symmetry in Fig.|21 left/right hand 
panel, respectively. However, the stellar dynamo-generated 
magnetic field outside the star develops in such a way that 
the azimuthal component of the field tends to be of one sign 
within each closed loop: within almost the whole 'dipolar' 
loop of Fig. (left hand panel), it is B^ > outside the 
star, and within the 'quadrupolar' loops of Fig. |2| (right hand 
panel), it is B^ < outside the star. This behavior of the 
toroidal field is possible because we do not impose any northern 
southern hemisphere symmetry for the rotation, but the angu- 
lar velocity evolves in a dynamical way within the star, at the 
stellar surface and outside the star. 

The 'oscillation' of the magnetosphere is not periodic and 
results in quick temporal changes of the magnetic field geom- 
etry with many different field configurations of mixed parity. 
Polarity reversals of the magnetic field occur mostly at time 
scales of less than a day. The irregular behavior of the mag- 
netic field geometry with time, including the magnetic polar- 
ity reversals, is likely to be due to the assumption of strong 
convection in our model of the star, which led us to choosing 
supercritical values of the dimensionless a dynamo number 
in the star: 1Z s ^ aT > 100 with an effective resistivity of a few 
times 10 -3 . Since we are in a highly nonlinear regime, the 
magnetic field can then be very irregular in the regime of an 



mass loss rate is M w - m 
the mass accretion rate varies over a large range, between 
M accr * lQ-^Mpyr- 1 and M accr ^lO^Moyr" 1 (see 
Appendix | Appendix B:| for the unit for M). 

In Fig. [3] a space-time diagram is shown for Model Ref. 
No clear latitudinal migration pattern of the radial magnetic 
field at the stellar surface can be seen. (The same is true for 
Model SmS, Model SDd and Model NDd.) 



3.3. Stellar dynamo versus anchored stellar dipole 

We recall that the star is a part of the computational domain. 
Therefore, in our stellar dynamo models no boundary con- 
ditions are imposed in or around the star other than the pre- 
scription of a mean-field dynamo in the star with a positive 
a effect in the upper hemisphere. This leads to a dynamical 
evolution of the magnetic field, angular velocity and poloidal 
velocity also within the star and at the stellar surface in the 
stellar dynamo models, whereby the stellar magnetic field 
strength and geometry are governed by the stellar dynamo 
parameters and Ma sta r. The resulting configuration of 
the dynamo-generated stellar magnetosphere in our reference 
model (Model Ref; cf. Fig. |2| and the top row of Fig. EH is 
highly time-dependent on rather short time scales of usually 
less than a day, including change of symmetry and polarity 
reversals, and varying magnetic field strength. (In order to 
achieve a stellar field strength as high as 1 kG and higher, 
one would need to increase the dynamo parameters such that 
the field oscillations would probably have unrealistically high 
frequencies.) Without restricting boundary conditions in the 
star, the rotation period of the star is rather low, having val- 
ues of a few days at the surface around the equator. As a con- 
sequence, an overall fast stellar wind develops in the stellar 
dynamo models. This wind is supersonic and super- Alfvenic, 
with terminal outflow speeds of a few hundreds of km/s. 
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Fig. 4. The absence of pronounced current sheets in a model with stellar dynamo (top row) in contrast to the case of an 
imposed stellar dipole field (bottom row). Colors/gray shades: temperature as in Figs [I] and black dashed line: poloidal 
Alfven surface; black arrows: poloidal magnetic field vectors (top row and bottom right; length of vectors is weighed with 
w) and poloidal velocity vectors (bottom left), respectively. Top: Model Ref at times t « 1171 days (left) and t « 1183 days 
(right). Bottom: Model S at time t « 158 days (cf. Figs 12 and 16 of vRB04.) Note that in Model S, collimation at the upper 
and lower boundaries (z > 0.8 and z < —0.8) is an artifact due to boundary conditions. 



In the magneto spheric models of vRB04, where a dipole 
magnetosphere is anchored in the star, the rotation profile in 
the star is fixed in time to rather large rotation periods; in 
Model S of vRB04, the minimum rotation period in the star 
is about 10 days. As a consequence, the whole stellar wind 
is slow in these magneto spheric models of vRB04, namely of 
the order of a few tens of km/s. Model S of the present paper 



differs in that the angular velocity as well as the azimuthal 
magnetic field can now freely evolve also within the anchor- 
ing region (i.e. within the star and around the stellar surface.) 
As a result, parts of the stellar surface have a rotation pe- 
riod of only a few days; they are located at mid-latitudes. 
These are the latitudes where a fast wind is launched from 
the star, with terminal speeds comparable to those of the fast 
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Fig. 5. Magnetic star-disc coupling and mass accretion onto the star in Model SDd with stellar and disc dynamos, showing 
two distinct states: when the star is magnetically connected to the disc (t w 464 days, left column) and when the star is 
magnetically disconnected from the disc it « 469 days, right column). The upper two panels show temperature color/gray 
scale coded, together with poloidal magnetic field lines and poloidal velocity vectors. The lower two panels show gas density 
color/gray scale coded, together with poloidal magnetic field lines and azimuthally integrated poloidal momentum density 
vectors. Note that in this stellar dynamo model, star-disc coupling is realized by the disc dynamo-generated magnetic field. 



inner disc wind (around 200 km/s). High latitudes rotate more 
slowly and the inner stellar wind from these latitudes remains 
slow (about 30 km/s); see the lower left panel of Fig.0] We 
will discuss the launching and acceleration mechanisms of 



the winds and the angular momentum transport in more de- 
tail in Sect.lT5l 

We checked that in Model S, the closed stellar field lines 
rotate nearly rigidly with the star. We also note that these field 
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Fig. 6. Radial profiles of the magnetic field strength (left hand panel) and the angular velocity (right hand panel) in the 
midplane at the times of Fig. (Model SDd), t « 464 days (solid lines) and t ~ 469 days (dashed lines). The vertical lines 
mark the stellar and inner disc radii. The solid curves correspond to a time when the field in most parts of the gap between 
the star and the disc is very weak, and hence disc matter - carrying Keplerian angular velocity from the disc throughout large 
parts of the gap - can accrete onto the star (connected state). The dashed curves correspond to the other extreme, a time when 
the field in the gap is strong, so that accretion and rotation in the gap are suppressed (disconnected state). The stellar surface 
rotation period changes from 4 to 10 days at the equator; the disc remains in Keplerian rotation. 



lines tend to show small temporal oscillations in the angular 
velocity. Whereas the angular velocity in the gap is around 
1 (so the corotating disc radius is also around 1), the outer- 
most closed field line that just penetrates the inner disc edge 
is rotating fastest with the Keplerian angular velocity of the 
inner disc edge (ft « 3); we return to this below in Sect. 13.71 
In contrast to the stellar dynamo models, the imposed dipole 
magnetosphere is oscillating close to periodically with a pe- 
riod of roughly 20 days. The oscillation period could be re- 
lated to the advection time through the disc, which is around 
20 days because the disc accretion velocity is around 10 km/s. 

Another important difference between the models with 
stellar dynamo and with anchored stellar dipole is the way 
in which the stellar field connects with the field in the disc. 
In both models, the fields from star and disc arrange them- 
selves in such a way that no X-point forms. In the anchored 
stellar dipole model this leads necessarily to the formation of 
pronounced hot current sheets above and below the midplane 
(vRB04). This is not the case in the stellar dynamo models, 
where the stellar field geometry is sufficiently more complex 
and able to adjust itself such that no current sheets form (cf. 
upper and lower rows of Fig.|4]>. 

The disc field and disc wind properties, like field strength 
and geometry, as well as wind speed and mass loss rate, are 
fairly similar. This is expected because the disc is modeled in 
the same way except for small changes in the disc dynamo 
parameters. 

3.4. Magnetic star-disc coupling and accretion 

In our stellar dynamo models, for the stellar magnetic field 
to be strong enough to penetrate the disc in the presence of a 
disc dynamo, one needs a very strong stellar dynamo to cre- 
ate field strengths of more than 1 kG. In this case, however, 
the stellar wind is too strong to allow for accretion, in par- 



ticular to allow for a funnel flow along the stellar field lines 
leading to polar accretion. Therefore, in our stellar dynamo 
models an efficient accretion flow is only possible in mainly 
equatorial direction. Also in the models with imposed dipole 
magnetosphere of vRB04, star-disc coupling by the stellar 
magnetosphere requires a stellar surface field that is stronger 
than 1 kG. But in contrast to the stellar dynamo models, most 
of the stellar wind remains weak so that polar accretion along 
the magnetospheric field lines is possible in the stellar dipole 
models. 

Even in the absence of a disc field, a stellar field strength 
of at least about 500 G is necessary in order to produce a cou- 
pling between star and disc by the stellar dynamo field; this 
is discussed further in Sect. 13.51 This leads to sufficient stel- 
lar wind speeds of a few hundred km/s, but suppresses the 
already weak accretion. Weaker stellar fields produce slower 
stellar winds. However, the magnetic coupling to the disc is 
then essentially absent. 

In order to get a coupling of star and disc, we made the 
disc dynamo stronger. In Model SDd, the disc dynamo pa- 
rameters are doubled compared to our reference model (Model Ref), 
leading to rms disc field strengths that are about 1.5 times 
larger (around 40 G). At the same time, the a effect in the 
star is reduced by a third, which leads to somewhat smaller 
rms stellar field strengths (< 200 G). 

Figure|5]shows two different states regarding the star-disc 
coupling. On the left hand side, star and disc are magnetically 
connected by field lines generated by the disc dynamo, which 
allows for equatorial accretion of disc material along those 
lines onto the star; these field lines are very weak in the gap 
between the disc and the star, where the plasma beta is above 
100. The top accretion rates are around M accr w (2 ... 6) x 
1O -T M0 yr _1 , which is comparable to the minimum of the 
total stellar and disc wind mass loss rates, the latter assum- 
ing values of around M w i n d « (6 . . . 18) x 1O _T M yr _1 . 
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Table 2. Results for accretion and wind properties of the main models discussed in this paper. Note that we are interested in 
matter that is accreted onto the star, so the mass accretion rate M accr is given near the stellar radius. Since accretion is highly 
episodic and irregular in time, the values of M a ccr are peak values. Typical values for the accretion flow velocity are given in 
the disc as well as in the gap. Wind speeds are terminal speeds (in our computational domain, excluding possible boundary 
effects). M w i n d = M ste iiarwind + Miiscwind is the total wind mass loss rate. The accretion rate M accr and the wind mass 
loss rates Miliar wind and Mdisc wind are in units of M yr _1 , whereas the accretion flow velocity u accr and the wind speeds 
^stellar wind and ^discwind are in units of km/s. The asterisk indicates that the values are for inner and outer stellar winds in 
Model S. 



Model 


-^accr 


disc 
"'accr 


"'accr 


-^stellar wind 


^stellar wind 


-^disc wind 


^disc wind 


Maccr /M w ind 


N 


2 x l(T y 


10 


< 1 


4 x l(T y 


10 


(1.5... 2) x 1(T Y 


350 


1 : 100 


S 


(1...2.5) x 1(T 8 


5... 15 


15 


2 x 10" 8 


30 / 200 * 


(2... 3) x 1(T 7 


250 


1 : 10 


Ref 


(0.5... 1) x 1(T 7 


10 


60 


(1...2) x 1(T 7 


400 


(1...2) x 10" 7 


250 


1 : 4 


SDd 


(2... 6) x 1(T 7 


10... 20 


75 . . . 150 


(2... 6) x 10" 7 


450 


(4... 12) x 1CT 7 


250 


1 : 3 


NDd 


negligible 


1...3 


negligible 


3 x 10" 7 


300 


2 x 1(T 7 


200 


negligible 



The accretion flow velocity in the gap is rather high, up to 
75 km/s and occasionally up to 150 km/s, but still consider- 
ably smaller than free-fall. 

On the right hand side of Fig. |5J the magnetic field in 
the gap - again originating from the disc dynamo - is almost 
vertical and strong enough to act as a shield for the star, not 
allowing disc matter to cross field lines. This causes star and 
disc to be disconnected in any way. In this state, disc matter 
leaving the inner edge is diverted into the inner disc wind, 
and there is no accretion. 

In Fig. [6] we show radial profiles of the magnetic field 
strength and the angular velocity in the midplane of Model SDd 
at the same two times as in Fig.|3 The solid (red) curves cor- 
respond to a state when the field in the gap between the star 
and the disc is weak and hence the Keplerian angular velocity 
from the disc can extend almost all the way to the star. The 
other extreme corresponds to the dashed (blue) curves where 
the field in the gap is strong and the angular velocity in the 
gap is suppressed. 

Table|2| summarizes the results for the accretion and wind 
properties of the main models discussed in this paper. The 
main points are the following: First, in our models with stel- 
lar as well as disc dynamo, especially in Model SDd, the ac- 
cretion flow velocity reaches - for the first time in our set of 
models - values that agree much better with observations of 
around 100 km/s or higher. The larger accretion flow veloc- 
ity might be due to the much larger plasma beta in the accre- 
tion region in the gap. Second, in these models, the mass ac- 
cretion rate becomes comparable to the total wind mass loss 
rate. Third, the stellar wind mass loss rate becomes compara- 
ble to the disc wind mass loss rate in the presence of a stellar 
dynamo. Forth, the accretion flow through the disc is clearly 
slower when there is no disc field. This shows that a mag- 
netic torque in the disc is important as a driver of angular 
momentum transport outward, for the accretion flow through 
the disc to be efficient. Indeed, in the absence of a disc dy- 
namo, the disc wind does not carry much angular momentum 
(see Sect. 13.51) . and therefore there is less accretion through 
the disc. Accretion flow in the gap and mass accretion rate 
onto the star are negligible in Model NDd, because the stellar 
wind is already too strong. 



3.5. Driving and structure of stellar and disc winds 

As we have seen in Fig.|2| in the case with a stellar dynamo 
there is a strong stellar wind. In order to clarify its origin 
we need to look at the force balance. For the magnetic ac- 
celeration, we consider the Lorentz force and the gas pres- 
sure force, both along the poloidal field whose unit vector is 
B po \ = £? po i/|£?poi|, i.e. we consider the forces 



(5) 



pol -B pol -(JxB), F£l=-B pol -Vp. 

In an axisymmetric system, if the magnetic field were purely 
poloidal, the Lorentz force would point in the poloidal direc- 
tion that is perpendicular to B po \, and therefore it would be 
Fp^ as ^ = 0. This can be seen explicitly by writing down the 



detailed expression for F^^ in axisymmeuy 

F (mag) = Bpol _ [(v x B ^ x 



pol 
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0. Here, is the unit vector in 



which vanishes for 
azimuthal direction. However, in a rotating magnetized star- 
disc system there is of course a strong toroidal component of 
the magnetic field in certain regions, leading to F^^ being 
different from zero in those regions. In Fig. we plot the ratio 

I^kJT^ l/l^poi I - Because of how this ratio is defined, regions 
where this ratio is large indicate regions where magnetic ac- 
celeration is strong. 

Concerning the disc wind, note the light shades especially 
above and below the disc in the region of field lines originat- 
ing from the inner disc edge, indicating that magnetic accel- 
eration is strong there. This "conical shell" corresponds also 
to the region of largest angular momentum transport away 
from the disc, so that the conical inner disc wind is magneto- 
centrifugally accelerated along the field lines. This conical 
structure requires a disc dynamo that is saturated so that the 
field above and below the disc is well-organized. The struc- 
tured disc wind due to the disc dynamo has been discussed 
in detail in von Rekowski et al. (2003); it is essentially in- 
dependent of the modeling of the stellar field. No magneto- 
centrifugal driving of the inner disc wind occurs when the 
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Fig. 7. The force ratio, | i 7 ^™ | / 1 | , is shown everywhere, including the disc, but not in the inner part of the star. Regions 
with largest Lorentz force along the poloidal field correspond to regions of highest angular momentum transport outward, and 
thus to regions of magneto-centrifugal acceleration. These regions are also those with lowest temperature and lowest density. 
Plasma beta is low and wind speeds are high in these regions, as well. White dash-dotted: fast magnetosonic surface. Left. 
Model Ref. Right: Model S, where the outer stellar wind is very similar to the inner disc wind. 




Fig. 8. Model NDd (no disc dynamo) at a time t w 246 days. Note that if there is no "disturbing" disc field, then star-disc 
coupling is possible for a stellar field strength of a few hundred gauss, but there is virtually no accretion of disc matter onto 
the star. Left: Cf. Fig.[l] Right: Color/gray scale representation of the force ratio (also shown in the disc but not in the inner 
part of the star), together with poloidal magnetic field lines (white/red [gray] means right/left-handed). The white dash-dotted 
line is the fast magnetosonic surface. As in Hayashi et al. (1996), hot plasmoids are ejected with a stellar wind speed of about 
300 km/s which is more than twice the Keplerian speed at the inner disc edge. 
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disc is virtually non-magnetized, i.e. the disc dynamo is switchedhavior, and configurations favorable for magneto-centrifugal 



off (aQ lsc = 0); see Fig.[SJ right hand panel. Here, the regions 
above and below the disc are predominantly black and corre- 
spond to the regions with basically no angular momentum 
transport, which is suggestive of mostly pressure driving of 
the whole weakly magnetized disc wind. 

The nature of the stellar wind depends strongly on the 
field configuration of the stellar magnetic field. This can be 
seen in Fig. [8] which shows a time snapshot of a model where 
there is only a stellar dynamo; the disc dynamo is turned 
off. Above and below the disc, there are narrow structures at 
about 45° angle where magneto-centrifugal driving is impor- 
tant. A similar structure is also present in the model with an 
anchored dipole; see the right hand panel of Fig. Here, the 
magnetic field can act as an efficient lever arm that is required 
for acceleration. All magneto-centrifugally accelerated winds 
are supersonic but the Alfven surface is further away. Angular 
momentum is transported along field lines (originating from 
the disc) that make an angle away from the vertical axis that 
exceeds 30° (Blandford & Payne 1982). 



In Fig.|9|we plot the inverse plasma beta, (3 



± mag/-fgas 

\B\ 2 /(2/j j qp), where the symbol P gas = p has been used. It 
is noticeable that in the stellar dipole model, the inner stellar 
wind has very low plasma beta but no poloidal component 
of the Lorentz force along the field lines. This is because the 
magnetic pressure due to the poloidal field exceeds the gas 
pressure. 

In summary, it turns out that in the stellar dynamo mod- 
els, the overall stellar wind is driven by magnetic and gas 
pressure. This wind is hot, dense, very fast and rotating with 
a period P rot of a few days. However, in places where field 
lines assume a configuration such that they can act as a lever 
arm, magneto-centrifugal acceleration becomes possible. The 
magnetic field geometry of the stellar field generated by the 
dynamo in the star is complex with an irregular temporal be- 
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Fig. 10. Model SDd. T mat , n et(0 and T mag (r) at fixed spher- 
ical radius r = 0.35 (close to the star) for various times. The 
mean amplitude of the net material torque is comparable to 
that of the magnetic torque, but both torques are highly time- 
dependent and change sign in an irregular fashion. Time is 
normalized, with 302 and 305 corresponding to the times of 
Fig. |5J 464 days and 469 days, respectively. 



acceleration are assumed locally for short times; the exact 
behavior is sensitive to the stellar dynamo parameters. In the 
dipole model, the stellar wind is clearly structured. Here the 
inner stellar wind is mostly driven by poloidal magnetic pres- 
sure. This wind is hot, dense, slow and rotating with P rot ~ 
10 days. The outer stellar wind is magneto-centrifugally ac- 
celerated and very similar to the inner disc wind (cooler, less 
dense, fast and rotating with P rot ~ 3 days). The magneto- 
centrifugally accelerated regions of the winds generally cor- 
respond to the coolest wind regions and have temperatures 
between 10 4 and 10 5 Kelvin. In the dipole model, the outer 
stellar wind and the inner disc wind are separated by a hot 
current sheet of about 10 6 K. 

As remarked before, the imposed stellar dipole field yields 
very little toroidal field in the plasma above the star, giving 
a poloidal component of the Lorentz force along the field 
lines only in a narrow strip. A stellar dynamo produces a 
more irregular field, and there can be a significant gradient 
in the toroidal magnetic field, giving rise to the possibility 
©f magneto-centrifugal driving of cooler and less dense wind 
structures. 

3.6. Stellar spin-up or spin-down? 

In order to address the question of spin-up versus spin-down 
of the star, we need to calculate the sum of the various torques 
exerted on the star, i.e. the net material torque (T matj net)> me 
magnetic torque (T mag ) and the viscous torque. The net ma- 
terial torque is the difference in angular momentum flux to- 
ward the star, carried by matter, when comparing the angular 
velocity of the accreting or outflowing matter, ft, and the ef- 
fective angular velocity of the stellar surface, £7* = L*/I*. 
Here, L* is the total stellar angular momentum, I* = M*m 2 
is the star's effective moment of inertia, and w 2 is the effec- 
tive lever arm that is, for reasons that become evident in a 
moment, defined as a weighed average of the form 

(w 2 ) = j> qu vo \w 2 -dS I ' j> qu vo \ • dS, (7) 

where the integrals are taken over the stellar surface. This 
defines ft*, which replaces the intuitive definition of ft* as 
the stellar surface rotation rate (see Sect. 3.8 of vRB04 and 
the appendix therein). We assume that the effective lever arm 
is constant in time, which means that during accretion or mass 
loss of the star matter is redistributed in such a way that w 2 is 
constant. This is done because in the model the stellar interior 
is insufficiently resolved and any variation of w 2 would be 
artificial. It is therefore important that our spin-up or spin- 
down calculations are not affected by this. 

Given the definition ft* = L*/I*, we can calculate its 
rate of change as 

^* = L* / 1* — L*I* j I 2 . (8) 

Since we have assumed w 2 = const, we have = M*vu 2 . 
To determine the sign of Cl* we consider the sign of 

= L* - L*M*/M*, (9) 

where L* is taken to be the sum of material and magnetic 
torques; the viscous torque is neglected because we found 
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Fig. 9. The inverse plasma beta is also shown in the disc but not in the inner part of the star. Note the region of low plasma beta 
also in the region of the inner stellar wind in Model S (right hand panel). There, one has only the poloidal magnetic pressure; 
the toroidal field and the angular velocity are very small. Blue & Orange dashed: Alfven & sonic surfaces, respectively. Left. 
Model Ref. Right: Model S, where the outer stellar wind is very similar to the inner disc wind. 
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Fig. 11. Model SDd. T mat , n et(^, ©) and T mag (r, 6) at fixed spherical radius r = 0.35 (close to the star) as functions of 
latitude at two different times, t « 464 days (left: connected state) and t « 469 days (right: disconnected state); cf. Fig|3 Net 
material torque and magnetic torque are comparable in amplitude. Left: At this time, the latitudinally integrated net material 
torque (T mat , n et(^ = 0.35)) is positive and the latitudinally integrated magnetic torque (T mag)net (r = 0.35)) is negative. 
Low-latitude accretion leads to a material spin-up. Higher-latitude stellar wind leads to a material spin-down but also to a 
material spin-up in places. Right: At this time, T mditynet (r = 0.35) < and T mag5net (r = 0.35) < 0. There is no accretion of 
disc matter, and the stellar wind generally leads to a material spin-down. 



that it is much smaller than the other torques. Using L* = This term can be subsumed into the material torque by defin- 
Q*M*ml and M* = — § gu po \ ■ dS, the second term on the ing a net material torque as (cf. vRB04) 
right hand side of Eq. ^ can be written as 

-^~mat,net (0 T \ ^mat,net 

(r,9)d6 (11) 

Jo 

with 

- L*M*/M* = -O^M* = j> Q*gu pol w 2 • dS. (10) T mat?net (r, 6) = -27ruj 3 gu r (Q - ft*)- (12) 
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Fig. 12. Model SDd. Radial profiles of the latitudinally integrated torques (HTl and JT3l at the times of Fig. [H] t « 464 days 
(left: connected state) and t « 469 days (right: disconnected state). 



Here we have used spherical polar coordinates, (r, 6, so 
w = rsinO and dS = 27rrtz7dB. The magnetic torque is 
defined correspondingly, 

T mag (r) =r f T mag (r,9)d6, (13) 
Jo 

where 

T mag (r, 6) = 27rw 2 B r B^/fi . (14) 

Both radial and latitudinal torque profiles will be discussed. 
The sum of both magnetic and net material torques gives 

— ^mat,net(^*) H~ ^mag(^*)j (15) 

and its sign determines whether the star spins up or down. 
Positive means stellar spin-up (and, since w\ is constant in 
time, gain in mean specific angular momentum of the star, 
Z* = mlQ*), whilst negative means stellar spin-down (and 
loss in mean specific angular momentum of the star). 

In the following we first turn to the dynamo models, in 
particular to Model SDd. In Fig.[l0|we show, at a fixed spher- 
ical radius r = 0.35 (close to the star), the latitudinally in- 
tegrated torques T mat , n et(V) and T mag (r) for various times. 
Both torques are highly time-dependent and change sign in an 
irregular fashion. It turns out that close to the star, the mean 
amplitude of net material torque is comparable to that of the 
magnetic torque. 

At this point we need to mention that a positive (negative) 
net material torque can reflect two situations: matter rotating 
faster (slower) than the star is flowing toward the star, or mat- 
ter flowing away from the star is leaving (taking) rotational 
energy. 

In Fig. [H] we plot latitudinal profiles of the torques both 
during the connected state (left hand panel) and during the 
disconnected state (right hand panel). Note that during the 
connected state, i.e. when matter from the disc falls onto the 
equatorial regions of the star, there is a significant positive 
net material torque in a thin equatorial strip. This means that 
accretion leads to a material spin-up of the star at those lat- 
itudes. The stellar wind from higher latitudes can spin down 
the star or spin it up, depending on the ability of carrying 
away angular momentum with it. 



The magnetic field around the star is controlled by both 
stellar and disc dynamos. It is complex enough so that the re- 
sulting magnetic torque changes sign in latitude several times 
(cf. Torkelsson 1998). At mid and high latitudes, the mag- 
netic torque exerted on the star is due to the field generated 
by the stellar dynamo; this stellar field is not connected to the 
disc but can be positive as well as negative. In the equato- 
rial regions, the disc dynamo-generated field connects the in- 
ner disc regions to the star; here, the magnetic torque exerted 
on the star is negative which means that the disc field spins 
the star down. However, this magnetic spin-down torque can 
barely counteract the material spin-up torque. This is true also 
for the latitudinally integrated torques, showing a total spin- 
up of the star in the connected state (of Fig. [12| left hand 
panel) due to accreting matter, although the magnetic torque 
is predominantly negative. (The stellar surface lies effectively 
around 0.33, due to a smoothed profile for the star.) 

During the disconnected state, the stellar wind in total car- 
ries away angular momentum which leads to a material spin- 
down of the star. In this state, spin-down is supported by a 
negative magnetic torque; however, the field is not connected 
to the disc at all. In Model SDd, the overall angular mo- 
mentum transport by the stellar wind decreases further away 
from the star, so that T mat , n et becomes positive while T mag 
changes sign twice. However, the stronger the stellar dynamo 
is, the more efficient is the braking by the stellar wind. 

In conclusion, stellar spin-up due to accretion cannot al- 
ways be compensated by the magnetic torque when the lat- 
ter is negative. Stellar spin-down by a stellar wind is equally 
important to stellar spin-down by a magnetic field. The to- 
tal (latitudinally integrated) torques at the stellar surface are 
changing sign over time scales ranging from a few days up 
to about 30 days, so that the star is alternately spun up and 
down. 

In Fig. [El we show the torques for the model with the an- 
chored stellar dipole field, Model S. Note that here, unlike the 
stellar dynamo case, the magnetic torque exceeds the net ma- 
terial torque and the stellar magnetic field leads to a spin-up 
of the star. Fig.[E]also shows that the fast outer stellar wind - 
emanating from fast-rotating latitudes with a rotation period 
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Fig. 13. Latitudinal (left column) and radial (right column) profiles of the torques for Model S at t « 158 days. For the 
latitudinal profiles we have plotted the torques (fT2l and (TUT) at r = 0.33, and for the radial profiles we have plotted the 
latitudinally integrated torques (fTTt and (IT3l . 
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Fig. 14. Top row: Model ND at t ~ 376 days. In this model, the computational domain extends to [—2, +2] in z. Top left: 
Black arrows: azimuthally integrated poloidal momentum density vectors (not shown in the star); colors/gray shades: gas 
density. Top right: Poloidal magnetic field strength I?poi as a function of spherical radius r at various latitudes (0, 30 and 60 
degrees), compared with fields decaying as B* (r/r*)~ 1,3 . Bottom row: Poloidal magnetic field strength £? po i as a function of 
cylindrical radius w at various heights (z = 0.15, z = 0.45 and z = 0.75), compared with fields decaying as tu -1,3 . Bottom 
left: Model Ref at t « 1171 days. Bottom right: Model S at t ~ 158 days. The vertical lines mark the stellar radius (top right 
and bottom row) and the inner disc radius (bottom row). 



of about 3 days - does yield a negative net material torque 
at these latitudes, causing a material spin-down of these lat- 
itudes. The major part of the stellar surface, however, is al- 
ready rotating rather slowly, with a rotation period of about 
10 days. Moreover, the magnetic torque is much larger than 
the net material torque and in total spins the star up. 



3.7. Collimation of stellar and disc winds 

A formal way of determining whether or not there is collima- 
tion was proposed by Fendt & Cemeljic (2002). Following 
their proposal, we choose a box within each of the two winds 
and define the mass load onto the stellar and disc winds as 

/\>star f*n'J ^TTW QU^ dz 
fstar _ 1V1 vj _ Jz ~v-b 
J load 



=0.7 



Mstar 



=0.7 
zu=0.5 



2hwqu z dm 



(16) 



=0.7 
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and 



^disc 
J load 



AC 



disc 



=0.75 
=0.6 



2-kwqu^ dz 



j; 



^7=1.05 



2nujgu z dm 



(17) 



=0.75 



Mass loads less than one indicate collimation while mass 
loads greater than one indicate decollimation. We have cal- 
culated these quantities for Models Ref and S; see also Figs|2| 
and |4| (bottom left), respectively. For Model Ref, the results 
are slightly different during the predominantly dipolar phase 
(Fig.0 left hand panel, where f**™ « 0.5 and « 2. ..3) 
and during the predominantly quadrupolar phase (Fig. right 
hand panel, where « 0.7 and « 6.5). For Model S, 
we have f^ r d « 0.6 and « 0.9 (see Fig. El bottom 
left). These results suggest that at low heights, the stellar wind 
shows a stronger tendency toward collimation than the disc 
wind. Collimation of the disc wind requires a well-organized 
and sufficiently strong magnetic field above and below the 
disc. 

The poloidal collimation properties of the stellar and disc 
winds depend on how rapidly the poloidal magnetic field de- 
cays with cylindrical radius (Spruit et al. 1997). The magnetic 
pressure from the poloidal field outside the jet is important. 
Spruit et al. (1997) derived the criterion that B po \ oc 
with /i < 1.3 is sufficient for poloidal collimation. 

For an ordinary dipole, the (poloidal) field decays with 
distance like 1/r 3 . This is also true of a dynamo-generated 
magnetic field in a vacuum, but here we are dealing with a 
dynamo in a conducting medium. In order to determine the 
radial decay of a dynamo-generated stellar field in our simu- 
lations, we consider Model ND with no disc; see Fig. [21 top 
row. The plot in the upper right panel shows the (spherically) 
radial decay of the poloidal field at different latitudes. There 
is no clean power law behavior but it is clear that the decay 
is slower than 1/r 3 . In fact, the decay of the poloidal field is 
similar to 1/ r 1,3 just outside the star. 

The tendency to poloidal collimation of the inner stellar 
wind in a stellar dynamo model can be seen in Fig. [21 bottom 
left, where we have plotted the poloidal magnetic field as a 
function of cylindrical radius at the disc surface and in the 
corona. The inner disc wind shows a tendency to poloidal 
collimation over a broader range in w only at larger z. 

In Model S, the decay of the magnetosphere is B po \(r) ex 
1/r 3 except for the equatorial region where it is compressed 
by the disc field and therefore decaying more slowly. The 
decay with cylindrical radius is B po \(zu) cx 1/vu 1 - 3 in both 
magneto-centrifugally accelerated winds; see Fig. [21 bottom 
right. 

A general feature of our models is that our dipolar disc 
dynamo-generated field advected into the corona satisfies the 
criterion for poloidal collimation of Spruit et al. (1997) in a 
certain range over w and z. Our models cannot be compared 
directly with M. v. Rekowski et al. (2000), whose disc dy- 
namo with dipolar polarity satisfies the condition everywhere 
at the disc surface, because in their model there are no rever- 
sals in B z in the disc. 

Figure[21 left column, shows that in a model with stellar 
dynamo, mid-latitudes rotate more slowly than the equator 



around the stellar surface. There is latitudinal as well as neg- 
ative radial rotational shear also outside the star, leading to a 
magnetic field that is decaying more slowly than a dipole with 
distance from the star (see Fig. [21 top right). A similar be- 
havior in the angular velocity is found in Model Ref (a model 
with stellar dynamo, disc and disc dynamo). In a model with 
anchored stellar magnetosphere (Fig. [21 right column), the 
radial shear is about around r anc hor at low latitudes but 
there is positive radial shear around r anc hor at higher latitudes 
where the fast magneto-centrifugally accelerated outer stellar 
wind is launched. Around r = r anc h or , these are also the lat- 
itudes that rotate fastest. This shows again that fast-rotating 
latitudes are more likely to produce a fast wind. 

4. Conclusions 

The present work has shown that a dynamo-generated stellar 
field, instead of just an anchored stellar dipole field, can ren- 
der the interaction between protostar and its surrounding disc 
in several ways more realistic. First, the observed structure 
of protostellar fields is known to be more complex than just 
a pure dipole. Second, dynamo-generated stellar fields can 
contribute quite efficiently to launching an overall fast stel- 
lar wind (up to 450km/s terminal velocity, which exceeds 
the terminal disc wind speed that itself is about twice the 
Keplerian speed of the inner disc edge). Such behavior was 
never found with just an imposed stellar dipole field. A stellar 
wind is equally important as the magnetic field for spinning 
down the star - especially because the braking via the mag- 
netic interaction with the surrounding disc is now found to be 
much less efficient than previously thought. We also suggest 
that this stellar wind could contribute to the observed jets, to- 
gether with the wind from the inner part of the surrounding 
accretion disc. 

As before in the case of an imposed central dipole field, 
the interaction between star and disc tends to occur in a time- 
dependent fashion where the system alternates between phases 
where star and disc are magnetically connected and mass can 
flow onto the star, and when they are magnetically discon- 
nected and mass transfer is blocked. Unsteady star-disc in- 
teraction is now becoming a prominent feature of magneto- 
spheric accretion onto T Tauri stars seen both in observations 
(Bouvier et al. 2003, 2004) as well as in numerical simula- 
tions (Goodson et al. 1997, 1999, Goodson & Winglee 1999, 
Matt et al. 2002, vRB04). At the same time, inefficient star- 
disc coupling is beginning to be a well recognized feature 
when stars couple to the differentially rotating disc (Matt & 
Pudritz 2004, 2005). 

However, time variability is no longer periodic but more 
irregular and on small scales, which is due to the more com- 
plex magnetic field configuration. The episodic magnetic star- 
disc coupling is no longer facilitated by the stellar magneto- 
sphere but by the disc dynamo-generated magnetic field so 
that the episodic accretion flow is not a slow polar funnel 
flow but happens at low latitudes along a very weak advected 
disc field, where it is much faster. Stellar braking is not only 
due to a magnetic field but also due to a fast overall stel- 
lar wind, which develops in the presence of a stellar dynamo 
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Fig. 15. Note the presence of radial and latitudinal shear in a model with stellar dynamo and no disc (Model ND; left column), 
which is compared to Model S (right column). (A behavior of Vt similar to that of Model ND is found in Model Ref.) Top 
row. Angular velocity as a function of spherical radius at various latitudes (0, 30 and 60 degrees), averaged over times 307 to 
461 days (left) and at time t « 158 days (right), respectively. Bottom row. Angular velocity as a function of latitude at various 
spherical radii (stellar surface, 2 times stellar surface and 3 times stellar surface), averaged over times 307 to 461 days (left) 
and at time t « 158 days (right), respectively. 



and is driven by gas and magnetic pressure, with magneto- 
centrifugal acceleration at times and in places where the stel- 
lar field geometry is appropriate. This suggestion is supported 
by recent observational evidence for high speed stellar winds 
from T Tauri stars (Dupree et al. 2005). 

Extension to nonaxisymmetry has been made for some 
dynamo models, and the results will be reported in a forth- 
coming paper. Possible steps to verify the dynamo models 
further could be to relax the constraints imposed by the adop- 
tion of mean-field theory and by the restriction to a piecewise 
poly tropic model. Another serious restriction that should be 
relaxed in future models is the assumption of time-independent 
profiles for the shapes of star and disc. Ideally, such profiles 
should be calculated in a self-consistent manner. In this pa- 
per, we have focused on the wind launching mechanisms and 
accretion processes. Studying the shaping and collimation of 
the winds further away from the star-disc system, requires a 
bigger computational domain. 
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Appendix A: The piecewise poly tropic model 

As described in our earlier papers, we model a cool, dense 
disc embedded in a hot, rarefied corona by prescribing con- 
trasts of the specific entropy s such that s is smaller within 
the disc and larger in the corona. For a perfect gas this im- 
plies p = Kg 1 (in dimensionless form), with the polytrope 
parameter K = e s l Cv being a function of position. The same 
idea is also applied to the star. In the models considered here, 
we have an intermediate value for the specific entropy within 
the star, resulting in a relatively cool and dense star. We pre- 
scribe the polytrope parameter to be unity in the corona and 
less than unity in the disc and in the star, so we put 

e 3/Cp = 1 - (1-AiiscKdisc - (l-&tar)£star, (Al) 

where £disc and £ sta r are explained in Sect. 2 of von Rekowski 
et al. (2003). The free parameters, < /3disc , Pst&r < 1, 
control the specific entropy contrasts between the disc and 
corona, and between the star and corona, respectively. With 
prescription JAIL (3di sc — 1 means that the disc is absent. 

In the absence of the stellar and disc dynamos, our initial 
state is a hydrostatic equilibrium with no poloidal velocity, 
where we assume an initially non-rotating, hot and rarefied 
corona that is supported by the pressure gradient. We solve 
the vertical and radial balance equations, as described in von 
Rekowski et al. (2003). Since we model a disc that is cool and 
dense, the disc is mainly centrifugally supported, and as a re- 
sult it is rotating slightly sub-Keplerian. The resulting stellar 
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rotation is differential, and the initial stellar surface rotation 
periods depend on the chosen stellar radius. 

The temperature ratio between disc and corona is roughly 
/^disc- Assuming pressure equilibrium between disc and corona, 
and p ex pT for a perfect gas, the corresponding density ra- 
tio is then f3^ c . Thus, the entropy contrast between disc and 
corona chosen here ((3^ sc = 0.005), leads to density and in- 
verse temperature ratios of 200 : 1 between disc and corona. 

A rough estimate for the initial toroidal velocity, u^o, in 
the midplane of the disc follows from the hydrostatic equi- 
librium as u^o w y/1 — /^disc^K- For (3^ sc = 0.005, the 
toroidal velocity is within 0.25% of the Keplerian speed. 

Appendix B: Physical units of our variables 

The resulting velocity unit is [u] = c s o = 10 2 kms _1 , the 
length unit is [r] = GM*/[u} 2 « 0.1 AU, the time unit is 
[t] = [r]/[u] w 1.536 days, the unit for the angular velocity 
is [Q] = [u]/[r] « 7.5 x 10~ 6 s _1 , the unit for the rota- 
tion period is [P ro t] ~ 10 days, the unit for the kinematic 
viscosity and magnetic diffusivity is [v] = [rj] = [u] [r] « 
1.5 x 10 19 cm 2 s _1 , the unit for specific entropy is [s] = 
c p = 7/(7 - ~ 3.5 x 10 8 cm 2 s- 2 K-\ the unit 

for specific enthalpy is [h] = [u] 2 = 10 4 km 2 s -2 , the tem- 
perature unit is [T] = [h]/[s] « 3 x 10 5 K, the density 
unit is [g] = So/M ~ 7.5 x 10 _13 gcm -3 , the mass unit 
is [M] = S [r] 2 « 2 x 10 24 g « 1O- 9 M , the pres- 
sure unit is [p] = (7 — l)/7[^][/i] « 30gcm -1 s -2 , the 
magnetic field unit is [B] = [n](47r[^]) 1 / 2 « 30 G (it is 
po = 47r), and the unit for the magnetic vector potential is 
[A] = [B][r] « 4 x 10 13 Gcm. Further, the unit for the 
mass accretion rate and mass loss rates due to the stellar and 
disc winds is [M] = E [tx][r] « 2 x 10" 7 M o yr~\ the unit 
for the torques is [T] = GM*£ MM « 2 x 10 38 erg, and 
the unit for the luminosity is [C] = GM*E [w] = [T]/[t] « 
1.5 x 10 33 ergs _1 . [From this follows that even for a very 
low mass accretion rate of M accr = 1O _1O M yr _1 and a 
protostellar radius of r* = 5Rq, one finds an accretion lu- 
minosity that is as high as C accr = GM*M accr /r* ^ 2.5 x 

10 30 ergg -lj 



References 

Agapitou, V., & Papaloizou, J. C. B. 2000, MNRAS, 317, 273 
Blandford, R. D., & Payne, D. R. 1982, MNRAS, 199, 883 
Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 
169 

Bouvier, J., Dougados, C, & Alencar, S. H. P. 2004, Ap&SS, 292, 
659 

Brandenburg, A., Nordlund, A, Stein, R. R, et al. 1995, ApJ, 446, 
741 

Cameron, A. C, & Campbell, C. G. 1993, A&A, 274, 309 
Doble r, W., Stix, M., & Brandenburg, A. 2005, ApJ (submitted), 

arXiv:astro-ph/041 0645 
Dupree, A. K., Brickhouse, N. S., Smith, G. H., et al. 2005, ApJ, 

625,L131 

Durney, B. R., De Young, D. S., & Roxburgh, I. W. 1993, Sol. Phys., 
145, 207 

Fendt, C, & Cemeljic, M. 2002, A&A, 395, 1045 



Gomez de Castro, A. I., & Verdugo, E. 2003, ApJ, 597, 443 
Goodson, A. P., & Winglee, R. M. 1999, ApJ, 524, 159 
Goodson, A. P., Bohm, K.-H., & Winglee, R. M. 1999, ApJ, 524, 
142 

Goodson, A. P., Winglee, R. M., & Bohm, K.-H. 1997, ApJ, 489, 
199 

Hartmann, L., Hewett, R., & Calvet, N. 1994, ApJ, 426, 669 
Hawley, S. L., Reid, I. N., & Gizis, J. 2000, in From Giant Planets to 

Cool Stars, ed. C. A. Griffith & M. S. Marley (ASP Conference 

Series, Vol. 212), p. 252 
Hayashi, M. R., Shibata, K., & Matsumoto, R. 1996, ApJ, 468, L37 
Hirose, S., Uchida, Y., Shibata, K., et al. 1997, Publ. Astron. Soc. 

Japan, 49, 193 
Konigl, A. 1991, ApJ, 370, L39 

Krause, R, & Radler, K.-H. 1980, Mean-Field Magnetohydrodyna- 

mics and Dynamo Theory, Akademie-Verlag, Berlin 
Johns-Krull, C. M., & Valenti, J. A. 2001, ApJ, 561, 1060 
Johns-Krull, C. M., Valenti, J. A., Hatzes, A. P., et al. 1999a, ApJ, 
510, L41 

Johns-Krull, C. M., Valenti, J. A., & Koresko, C. 1999b, ApJ, 516, 
900 

Lovelace, R. V E., Romano va, M. M., & Bisnovatyi-Kogan, G. S. 

1995, ApJ, 275, 244 
Matt, S., Goodson, A. P., Winglee, R. M., et al. 2002, ApJ, 574, 232 
Matt, S., & Pudritz, R. E. 2004, ApJ, 607, L43 
Matt, S., & Pudritz, R. E. 2005, MNRAS, 356, 167 
Meinel, R., & Brandenburg, A. 1990, A&A, 238, 369 
Menard, R, & Duchene, G. 2004, A&A, 425, 973 
Miller, K. A., & Stone, J. M. 1997, ApJ, 489, 890 
Nordlund, A., Brandenburg, A., Jennings, R. L., et al. 1992, ApJ, 

392, 647 

Papaloizou, J. C. B., & Terquem, C. 1999, ApJ, 521, 823 
Pudritz, R. E. 2004, Ap&SS, 292, 471 

von Rekowski, B., Brandenburg, A., Dobler, W, et al. 2003, A&A, 
398, 825 

von Rekowski, B., & Brandenburg, A. 2004, A&A, 420, 17 (vRB04) 
v. Rekowski, M., Rudiger, G., & Elstner, D. 2000, A&A, 353, 813 
Romanova, M. M., Ustyugova, G. V, Koldoba, A. V, et al. 2004, 

ApJ, 616.L151 
Shu, R, Najita, J., Ostriker, E., et al. 1994, ApJ, 429, 781 
Spruit, H. C, Foglizzo, T., & Stehle, R. 1997, MNRAS, 288, 333 
Stempels, H. C, & Piskunov, N. 2002, A&A, 391, 595 
Tavakol, R., Tworkowski, A. S., Brandenburg, A., et al. 1995, A&A, 

296, 269 

Tobias, S. M., Brummell, N. H., Clune, T. L., et al. 1998, ApJ, 502, 
L177 

Torkelsson, U. 1998, MNRAS, 298, L55 
Urpin, V, & Brandenburg, A. 1998, MNRAS, 294, 399 
Uzdensky, D. A., Konigl, A., & Litwin, C. 2002, ApJ, 565, 1191 
Valenti, J. A., & Johns-Krull, C. M. 2004, Ap&SS, 292, 619 



Id: paper.tex,v 1.120 2005/09/10 12:23:03 gitta Exp 



